Griffiths phases in the contact process on complex networks 
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Dynamical processes occurring on top of complex networks have become an exciting area of re- 
search. Quenched disorder plays a relevant role in general dynamical processes and phase transitions, 
but the effect of topological quenched disorder on the dynamics of complex networks has not been 
systematically studied so far. Here, we provide heuristic and numerical analyses of the contact 
process defined on some complex networks with topological disorder. We report on Griffiths phases 
and other rare region effects, leading rather generically to anomalously slow relaxation in general- 
ized small-world networks. In particular, it is illustrated that Griffiths phases can emerge as the 
consequence of pure topological heterogeneity if the topological dimension of the network is finite. 

PACS numbers: 



I. INTRODUCTION 



The study of complex networks has been nourishing since the introduction of a simple model describing the emer- 
gence of scaling in random networks Multidisciplinary applications involve, for example, the World Wide Web, 
various biological, sociological and technological networks. These are, in many cases, scale-free networks (for recent 
reviews see Other families of complex network models are those composed of a d-dimensional regular lattice 

and additional long edges These arise e.g. in sociophysics in the context of conductive properties of linear 
polymers with cross-links that connect remote monomers [f| , in public traffic systems @ , in the case of nanowires Q , 
and in many other examples. In general, a pair of nodes separated by the distance I is connected by an edge with the 
asymptotic probability for large I: 

p{i) = pr s . (i) 

In the special case s = 0, the edges exist with a length-independent probability, as in small world networks, therefore 
we call the nets with generic s generalized small world networks (GSW). If s > 2, they are characterized by a finite 
topological dimension D, i.e. N(l) ~ l D , where N(l) is the number of nodes within a distance I from a given node. 
(Note that " small- world" property means that the number of neighbors of a node increases exponentially with the 
distance from it, i.e. formally D = oo). 

Dynamical processes defined on networks are of recent interest, for example in the context of optimization dynamics 
of spreading or transport processes @- Different dynamical processes, defined on regular lattices, often exhibit scaling 
behavior, that can be classified into basic universality classes @. Instead, the existence of scaling universality classes 
on network architectures is not clearly established. In many cases, the presence of short average distances induces 
mean-field type of (fast) dynamics. It is also well known that structural heterogeneity may lead to more complex 
behavior; for example, on scale- free networks, characterized by a power-law degree distribution P(k) oc fc -7 
critical exponents may vary with the degree exponent 7 [1, [H|. In related 'annealed networks', where the links change 
rapidly, 7-dependent critical exponents have also been reported [ToL [111]. Real networks are, however, 'quenched' in 
many cases, i.e. the topology changes slowly with respect to the dynamical process evolving on them. 

In the context of the statistical physics of models defined on Euclidean lattices, it is well known that quenched 
randomness can generate "rare- region effects" in the so called Griffiths phase (GP), where, in many cases, algebraic 
scaling is observed (as opposed to pure systems in which power-laws are observed only at critical points) with scaling 
exponents changing continuously with the control parameter (l3l - fl6l |30| . Furthermore, at the phase transition the 
evolution becomes logarithmically slow, and it is controlled by 'activated scaling'. A nontrivial question can be posed: 
Under which conditions do rare-region effects and GPs occur in network systems? Can they emerge out of topological 
disorder alone (i.e. without disorder in the transition rates)? These problems have not been systematically studied 
so far. At first sight one might guess that in network models the dynamics should be very fast, due to the strong 
inter-connectedness, i.e. information can propagate very efficiently throughout the system, resulting in mean-field 
like, exponential relaxation. Recently, however, generic slow dynamics has been reported to occur in various network 
models; its origin has been attributed to heterogeneity [ItJ or to local bursty activity patterns [H, [l9[. To tackle 
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this type of problems, here we present a study of the contact process (CP) [20J (i.e. the simplest possible model 
for epidemic spreading, or information flow) on different GSW networks, including regular networks and non-regular 
ones. 

II. NUMERICAL ANALYSES 

We study the CP on random networks composed of a c?-dimensional lattice and a set of long edges with unbounded 
length. Any pai r of nodes, separated by the distance I, is connected by an edge with some ^-dependent probability, 
Eq.([T]) P. |24l. |25| . An intriguing feature of these type of graphs for d = 1 is that in the marginal case (s = 2) intrinsic 
properties show power-law behavior and the corresponding exponents vary continuousl y w ith the prefactor /3. Indeed, 
the topological dimension D of such networks has been conjectured to depend on (3 [2J]. Instead, for larger values 
of s long-ranged links are irrelevant, and D((3) — 1, while for s < 2 the topological dimension diverges. It has been 
claimed in a recent paper [28[ that if D((3) is finite, Griffiths phases and similar rare-region effects can appear. 

We have considered the CP [2fj], in which each infected (active) sites is healed at rate 1, whereas each of its 
nearest-neighbor sites is infected at rate X/k (where k is the degree of the site) so that the total infection rate is A. 
In numerical simulations we update active sites randomly and increment the time by 1/N a , where iV a is the number 
of active sites. In this way, all active sites are updated on average once every time unit. For a critical infection rate 
A c we find a phase transition from the active to an absorbing phase (2lj . with vanishing density of infected sites. 

Wc have performed numerical simulations initiated either from a fully active state or from a single active seed. In 
the former case, the density p(t) of active sites has been measured. 

A. Non- regular random networks 

The precise definition of the networks we consider is as follows [2^, [24|. We start with TV nodes, numbered as 
1,2, ... ,N and define the distance between node i and j as I — min(|i — j\, N — \i — j\). We connect any pair of 
sites separated by a distance I = 1 (i.e. neighboring sites on the ring) with probability 1, and pairs with I > 1 with 
a probability P(l) = 1 — exp(— j3l~ s ). This implies P(l) = j3l~ s at large distances. In the case s > 2, long-range 
links do not change the topological dimension, which remains unity, and hence the critical behavior is expected to be 
similar to that of the one-dimensional CP with disordered transition rates. In the latter model, the critical dynamics is 
logarithmically slow, for instance, starting from a fully active state, the density decays asymptotically as p cx (ln(i))~ Q 
with a = 0.38197 [lg. For s = 3 (and /3 = 2) wc have run extremely long (t max = 2 28 MCs) simulations on networks 
with number of nodes N = 10 5 . Besides observing a GP, we have found that data are compatible with the above 
form of logarithmic dependence on time, and have located the critical point at A c = 2.783(1). As Fig. QJa) shows 
the assumption on activated scaling with 5 = 0.38197 is satisfactory, although only after an extremely long crossover 
time to be discussed later (29j . 




(a) (b) 

FIG. 1: Time-dependence of the density in the s — 3, /3 = 2 random network. Inset: local slopes of the same (lower curves), and 
local slopes of the survival probability in the Id QCP simulations (dashed curves). Slow convergence of the effective exponents 
can be observed, (b) Phase diagram of QCP on the ER graph for r = 0. 



For the marginal case s = 2, where the topological dimension is a continuous function of j3, numerical simulations 
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indicate qualitatively different scenarios of the phase transition for different values of j3. If j3 is small enough, GP is 
observed and the scaling at the transition is of logarithmic form. However, when j3 is large enough, the GP seems 
to be lacking and the critical dynamics follows a power- law. The latter type of behavior is observed also for s < 2, 
where formally D = oo. 

An "annealed" counterpart of CP on the above random networks is the Id CP with Levy flight distributed activation 
probabilities P(r) oc r~ d ~ a . In this model the dynamical exponents are known to depend continuously on a [23]. 
The estimated dynamical exponent of the CP on the above random networks in case of power-law critical behavior is 
found to be compatible with the dynamical exponent of the Levy-flight CP with d=l. For further details the reader 
is referred to [28[ [3 . 

B. 3-regular (or cubic) random networks 

In the networks studied so far the degree of nodes was random. In the following we consider networks with a 
"weaker" topological disorder, in the sense that the degree of nodes is constant (3). Such random networks with 
nodes of degree 3 can be constructed in the following way [25|. A one-dimensional periodic lattice with N sites is 
given, where the degree of all sites is initially 2. Sites of degree 2 are called "free sites". Let us assume that N is 
even and k is a fixed positive integer. A pair of free sites is selected such that the number of free sites between them 
is fc — 1 (the number of non-free nodes can take any value) and this pair is then connected by a link. That means, 
for k = 1, neighboring free sites are connected, for k = 2 next-to-neighboring ones, etc. This step, which raises to 3 
the degree of two free sites, is then iterated until 2(k — 1) free sites are left. These are then paired in an arbitrary 
way, which does not affect the properties of the network in the limit N — > oo. In the resulting network, all sites are 
of degree 3, and one can show that the probability of edges is given by Eq. ([T]) with s = 2. For certain networks of 
this type it has been demonstrated that the long-ranged connections result in a finite topological dimension which is 
less than one (25|. Furthermore, in the case of aperiodic networks, the critical exponents of the CP depend on the 
underlying aperiodic modulation [26| . In the following, we concentrate on random networks in which the pairs to be 
connected are selected randomly (with equal probability) in the above procedure for a fixed k. In this case, one can 
show that the prefactor in Eq. (QJ is given by j3 = k/2. 

We have performed numerical simulations for the CP on k = 1 random networks with N = 10 7 nodes. The 
averaging was done over 200 different network realizations for each A value and the maximum time is t max < 2 26 . 
As Fig. HJa) shows, power-law decay with continuously changing exponent emerges for a range of A-s. By analyzing 
the effective decay exponents, defined as the local slopes of the density a c g(t) = — (In p(t) — lnp(i'))/(ln(t) — ln(i')), 
where t/t 1 = 2, the curves do not level-off for large times, instead one can observe a small drift. One can easily derive 
that the functional dependence of the effective exponent a c g(t) = a + x/ln(l/t) is related to the scaling correction 
p(t) = t~ a \n x (l/t). By plotting a c ff (<) on logarithmic time scales (see inset of Fig. ^a)) it turns out that the drift 
describes logarithmic corrections to power-laws and that a is a non-universal quantity which depends on A. The 
possibility of logarithmic corrections has already been pointed out in the case of CP with quenched disorder fTo| . 




(a) (b) 

FIG. 2: Density decay in 3-regular, random networks, (a) GP in k = 1 between A = 2.535 and A = 2.57 (from bottom to top), 
(b) Critical point for k = 2 between A = 2.156 and A = 2.1598 (from bottom to top). The inset shows the corresponding local 
slopes. 

In accordance with results on non- regular networks with s = 2, the phase transition is qualitatively different for 
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k = 2 which corresponds to a higher value of (3. Here, a standard phase transition seems to appear at A c = 2.15835(5), 
with the density decay exponent a ~ 0.51(5) (Fig. [UJb)). This value agrees again with that of the Levy flight CP 
with d = 1 and d + a = 2. 



III. QCP ON ERDOS-RENYI NETWORKS 



Apart from studying the effect of heterogeneous (i.e. disordered) topologies on the dynamics of the contact process, 
we have also scrutinized the behavior of the disordered contact process on random networks. In particular, we have 
considered CP on ER graphs [22j with a quenched disordered infection rate: a fraction 1 — q of the nodes (type-I) 
take value A and the remaining fraction q (type-II nodes) take a reduced value Ar, with < r < 1. Pair mean-field 
approximations [28l . I29I ] lead to the following critical threshold 

K(q) = TTT-^-r ~, ~ • (2) 

(k) - 1 1 - q 

Type-I nodes experience a percolation transition, where the type I-to- type I average degree is 1, i.e. at q per c = 1 — {fc}^ 1 - 
For q > q perc activity cannot be sustained: type-I clusters are finite and type-II ones do not propagate activity. 

Numerical simulations and optimal fluctuation theory for (k) = 3 (q per c = 2/3) leads to the complex phase diagram 
shown on Fig. [T^b). In agreement with ^ one finds a critical active/ absorbing phase transition, but below the 
percolation threshold the absorbing phase splits into parts. In particular, for A > 4.5 (iv) rare, percolating regions 
may occur, which cause power-law dynamics (i), i.e. a Griffiths phase, compared to the inactive phase of simple CP, 
which is purely exponential. Further details can be found in [2^, [29[ . 



IV. CONCLUSIONS 



We have illustrated that, besides quenched site disorder, topological disorder by itself can result in slow dynamics 
and non-universality, at least for the contact process. We expect these results to have a broad spectrum of implications 
for propagation phenomena and other dynamical process taking place on networks, and to be relevant for the analysis 
of both models and empirical data [13] ■ We have claimed in a different publication [28| that having finite topological 
dimension is a necessary condition for the occurrence of slow, GP type of evolution in complex networks, at least for 
CP type of dynamics. Investigation of other factors promises to be an interesting, open field of research. 
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